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The advantage of high efficiency, low SFC (Specific Fuel Consumption), ultra-high bypass ratio turbofan engine 
attracts more and more attention in modern commercial engine. The intermediate turbine duct (ITD), which con- 
nects high pressure turbine (HPT) with low pressure turbine (LPT), has a critical impact on the overall perfor- 
mance of turbine by guiding flow coming from HPT to LPT without too much loss. Therefore, it becomes more 
and more urgent to master the technique of designing aggressive, even super-aggressive ITD. 

Much more concerns have been concentrated on the duct. However, in order to further improve turbine, LPT 
nozzle is arranged into ITD to shorten low pressure axle. With such design concept, it is obvious that LPT nozzle 
flow field is easily influenced by upstream duct structure, but receives much less interests on the contrary. In this 
paper, numerical method is used to investigate the effects of length of ITD with upstream swirl blades on LPT 
nozzle. Nine models with the same swirl and nozzle blades, while the length of ITD is the only parameter to be 
changed, will be discussed. Finally, several conclusions and advices for designers are summarized. 

After changing axial length of ducts, inlet and outlet flow field of nozzle differs, correspondingly. On the other 
hand, the shearing stress on nozzle blade (suction and pressure) surface presents individual feature under various 


inlet flow. In addition to that, “Clocking-like effect” is described in this paper, which will contribute much to the 


pressure loss and should be paid enough attention. 
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Introduction 


It is pursued all the time to design an engine with 
higher thrust-weight ratio and lower cost, which requires 
ultra-high bypass ratio to lead larger diameter fan. How- 
ever, the tip speed of fan is limited to subsonic, so as to 
keep acceptable loss and be sure security. As the result, 
the rotation speed has to be decreased, so does the LPT 
rotors. With unchanged mass flow, the medium diameter 
of LPT will increase. So it becomes more and more im- 
portant, but more and more difficult to guide flow leav- 
ing HPT to LPT without severe pressure loss. The ques- 
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tion will be how to achieve it? ITD is believed to be in- 
vestigated essentially in more details, so as it could be 
designed more aggressively. More and more attention has 
been paid to this area in recent years. 

Before designing aggressive, even super- aggressive 
ITD, the primary work is to understand flow mechanism 
of ITD. Many researchers made hard efforts to investi- 
gate this duct, from the original simple annular diffuser 
to S-shaped duct, and modern aggressive duct finally. 

The investigation about flow in ITD began from sim- 
ple annular diffuser, which was done by Stevens et al.!"! 
at the University of Durham, and then simple S-shaped 
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Nomenclature and Abbreviation 


HPT High Pressure Turbine 

ITD Intermediate Turbine Duct 
LPT Low Pressure Turbine 

Ma Mach Number 

SFC Specific Fuel Consumption 


duct was measured and discussed by many other re- 
searchers'”!"“!, But subject to the limitation of measure- 
ment technology, the flow field gotten was crude and 
inaccurate. 

Miller et al.'°! described both the migration and dissi- 
pation of flow in the ITD with upstream high-pressure 
turbine stage in a transonic facility. The author found that 
the inlet flow to a downstream vane was dominated by 
two co-rotating vortices, the first caused by the rotor 
tip-leakage flow and the second by the rotor hub second- 
ary flow, which was the first to describe the detailed vor- 
tices development within ITD as the author knew. And 
then downstream nozzle vanes were presented in the duct 
to get a more realistic flow field. In a large-scale, 
low-speed, experimental turbine facility, Axelsson et 
al." investigated the flow development for different 
turbine conditions. 

In addition, CFD simulation, as an effective means, 
was adopted to investigate the ducts. Wallin et al.!*! and 
Osso et al.”! studied the flow field in an aggressive ITD 
under design inlet condition and off-design inlet condi- 
tion. They found that the total pressure loss increased 
under off-design condition, especially in the ducts with 
separation. Peter et al.''°! made experimental and numer- 
ical investigations of an aggressive intermediate duct 
configuration. With upstream HPT and downstream LPT 
nozzle, the author found that the hub and shroud pressure 
profiles were well predicted while the total pressure loss 
was over predicted in the numerical simulations. 

Steady flow field is important to understanding the 
mechanism of ITD, but whether the unsteady flow field 
is similar to the steady condition? With numerical simu- 
lation, Davis et al.''"! presented the unsteady flow in a 
transonic turbine stage, and compared their result with 
experiment, which showed overall good agreement be- 
tween the simulation and experiment data. 

Gottlich et al."*! built a transonic test turbine facility 
to explore flow within aggressive ITDs in Graz Univer- 
sity of Technology. With this facility, Gottlich et al.” 
studied the steady and unsteady flow field in an aggres- 
sive ITD downstream of a single-stage turbine. On this 
basis, a much shorter so-called super-aggressive interme- 
diate duct (20% shorter than the aggressive duct) was 
designed by Géttlich et al. ''*!, And Marn et al.l" dis- 
cussed the flow field in this super-aggressive ITD and 
compared it with aggressive ITD. With numerical simu- 
lation and experiment, the authors gained a more detailed 
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SP Static Pressure 

Cpt Total Pressure Loss Coefficient 
Pt Local Total Pressure 

Pti Inlet Total Pressure 


insight into the flow behavior of a fully separated duct. 

Except for investigating the steady and unsteady flow 
field within different ITDs, the influence of nearby tur- 
bine stages on ITD was researched. 

Dominy et al. discussed the influence of swirl!'®! and 
upstream wake!’”! on ducts. And then the effect of blade 
tip clearance!'®!!*7) inlet swirl angle distribution PP, 
upstream flow field’?! were reported by many research- 
ers. On the other hand, the flow fields in ITDs with dif- 
ferent geometries, such as various streamline curva- 
tures?°1?7) area ratios!?*!", diffusion rates°!, and so on, 
were described as well. Through different methods, in- 
vestigators got many meaningful conclusions, and gave 
many valuable advices for designing ITDs. 

According to the above introduction, muck work was 
finished. Steady, as well as unsteady flow field in annular 
diffuser, S-shaped ducts, and aggressive ITDs were de- 
scribed with experimental and numerical methods. The 
systematic knowledge about ITD flow field has been 
built. Besides flow field analysis in ITDs, the influences 
of blade tip clearing, upstream wake, geometry of area 
ratio, diffuser ratio, curvature, and so on, on the flow 
field of ITD were investigated. What’s more, through 
comparison with experimental results, a common conclu- 
sion was drawn that numerical simulation could get a 
good enough flow field and ITD performance, which 
could be used to analyze and evaluate the flow in ITD. 

Many excellent studies have been done about perfor- 
mance of ITD under different conditions, which are im- 
portant for designing aggressive diffuser duct. However, 
how nearby turbines are influenced by duct remains un- 
known, which is important for designing future turbine as 
well. This issue should be analyzed and discussed to 
draw clear conclusions and to give advisable guidance 
for designers essentially. 

In this paper, ITDs with the same upstream swirl 
blades and downstream LPT nozzle are designed, as 
shown in Figure 1. Nine models with different duct 
length are simulated with numerical method to investi- 
gate the effect of ITD length on LPT nozzle. As shown in 
Figure 2, four typical models are presented. The author 
hopes this paper could give an insight into how to design 
an ITD, as well as nearby turbine stages with more prop- 
er length while ITD is working under engine environ- 
ments. Meanwhile, findings revealed by the current paper 
will provide meaningful insights for further optimization 
of the ITD through the engineering perspective. 
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Fig. 2 The sketch of different length ITDs 


Numerical Model and Boundary Conditions 


The original aggressive ITD prototype in this paper 
comes from another investigation paper, which focused 
on the effect of rising angle on ITD performance. In that 
paper, eight ITD models with the same upstream swirl 
blades were simulated. Area ratio (=1.57) and non- 
dimensional length (L/H = 3.4) of ducts remained un- 
changed, while rising angle (8°~45°) was the only para- 
meter modified. On the basis of that paper, the ITD with 
25° rising angle is chosen to be the target object further 
investigated. 

After optimizing LPT nozzle blades to improve flow 
field again, the prototype model (Aero ratio = 1.57, 0 =2 
5°, L/H = 3.4) is defined. In order to understand how ITD 
length influences downstream LPT nozzle, nine models 
modified of L_3.4 (L/H = 3.4), L_3.2 (L/H = 3.2), L_3.0 
(L/H = 3.0), L_2.8 (L/H = 2.8), L_ 2.6 (L/H = 2.6), L_2.4 
(L/H = 2.4), L_2.2 (L/H = 2.2), L_2 (L/H = 2), L_1.8 
(L/H = 1.8) are simulated in parallel. Swirl blades and 
nozzle blades remain unchanged, as Table | lists. Blade 
numbers of swirl and nozzle are set the same 44 to guar- 
antee the results being credible, from which wake of 
swirl blade could be transported to LPT nozzle accurately 
without circumferential average. 


Table 1 Blade Parameters 


Swirl Nozzle 
Blade Number 44 44 

Axial Chord 21.64 28.88 
Exit Angle -20.80 61.74 
Inlet Angle 0 -31.39 
Stagger Angle -12.79 34.23 
Unguided Turning 0 15.46 

Pitch/Chord 0.77 0.77 

Tmax/Axial Chord 0.07 0.21 
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The boundary conditions for all models are the same, 
with a constant in-mass-flow of 8.52 kg/s at inlet, and 
back pressure of 240 kPa. So the Reynolds number is 
kept 3 x 10° to remove its influence on boundary layers. 

The commercial CFD software of Numeca is used. 
The topology of building model mesh is O4H, and 
matching periodicity technology is used to escape nu- 
merical dissipation so as to get an accurate flow field in 
blade passages. Figure 3 is the numerical model and 
structured grids of L_2.6 as an example created with IGG 
The mesh at blade root and leading edge is magnified 
locally to show the mesh detail. The orthogonality of 
mesh is good enough to reflect the flow field. And Y* 
near wall is set below | to fulfill the requirement of the 
Shear Stress Transport (SST) turbulence model. With the 
same count of blades, full non-matching frozen rotor 
between swirl blades and LPT nozzle is chosen to pro- 
vide steady upstream wake for ITD. 

The quantity of model mesh is an important factor for 
numerical results. So a test of mesh independence to si- 
mulation result is necessary. Figure 4 is the relationship 
between total pressure ratio (of model inlet to model out- 
let) and mesh quantity of L_3.4. This regularity applies to 
other models as well. From this chart, the total pressure 
ratio will not be influenced any more when the mesh is 
more than 3 million, which means that the flow field is 
independent of mesh count. So each model will be simu- 
lated with more than 3 million meshes. 


Fig. 3 Numerical Model and Mesh of L_2.6 
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Fig. 4 Test of Mesh Independence 
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Results and Discussion 


After simulating all the models, it is found that separa- 
tion occurs in ITD models shorter than L_2.2, especially 
L_1.8, in which separation is very severe at leading edge 
of LPT nozzle blade. Because all the nozzle blades re- 
main the same with L_3.4 so as to investigate the influ- 
ence of ITD length separately, blade profile is not opti- 
mizational except L_3.4. 

In the following discussion, four typical models of 
L 3.4, L 3, L 2.6, L_2.2 are chosen to be analyzed to 
draw objective conclusions. As what has been stated 
above, this paper focuses on the effects of ITD length on 
LPT nozzle, so flow field and performance of nozzle 
blades will be discussed in detail, while flow field within 
ITD is paid less attention. Inlet as well as outlet angle 
distribution will be discussed preferentially. 

The radial yaw angle distribution at LPT nozzle inlet 
differs when the length of upstream ITD changes. As 
shown in Figure 5, the overall trends for all the ducts are 
very similar. The peak magnitude locates at about 95% 
passage height, and two relative local troughs appear at 
about 19% and 78% height, which mean that the ITD 
length will not change the radial angle distribution trend 
too much. 


LP Inlet Angle 


15 20 25 30 35 40 45 
Yaw Angle 


Fig.5 Radial yaw angle at LPT nozzle inlet 


But local detail is various. An evident phenomenon is 
that the regions near walls are more easily influenced. 
According to Figure 5, the deviation amplitude difference 
near walls is large. A circumferential flow about 42 de- 
grees in model L-2.2 compared with 32 degrees in L-3.4 
presents a severer deviation near outer wall. Tangential 
velocity is larger than axial velocity, which means thicker 
boundary layer on outer casing in shorter models. On the 
other hand, within the main flow near mid-height passage, 
the angle is more stable relatively with a more axial flow 
in shorter ITDs. And this trend is severer near hub. So the 
boundary layer on hub of shorter duct is thinner corres- 
pondingly. Obviously, the attack angle of LP nozzle blade 
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will be more negative in more aggressive ITD model, 
which may induce separation on the pressure surface to 
increase pressure loss. This phenomenon will be de- 
scribed in following contents. Though the boundary layer 
distribution within ITD is not emphasized in this paper, 
such conclusion should be acknowledged that boundary 
layer on hub will be thinner, while it is thicker on casing 
when ITD length decreases. Without doubt, the blade 
profile of LPT nozzle near walls should receive more 
attention in order to control pressure loss. 

After developing through the same nozzle passage, 
flow information at the outlet of LPT nozzle shows the 
effect of upstream ITD is opposite, as Figure 6 shows. 


LP Outlet Angle 


Passage Height 


Yaw Angle 


Fig.6 Radial yaw angle at LPT nozzle outlet 


Near walls, the flow yaw angles among various ducts 
have similar trend relatively that the number and located 
height of local extremums are almost the same, which is 
identical with the distribution at inlet. The angles near 
walls are larger within boundary layers, while more axial 
in the main stream, especially at about 15% and 75% 
passage height where upper passage vortex and lower 
passage vortex locate. 

However, the radial angle details show individual fea- 
tures above 15% to 80% passage height. The regularity 
where yaw angle changes within different ducts shows 
that the peak magnitude migrates with a period radially 
in Figure 6. Actually, the angle distribution is influenced 
mainly by passage vortices and flow near blade corners. 
With the same blade profile and downstream passage, the 
pressure field near hub remains almost unchanged, given 
the thinner boundary layers at nozzle inlet. While near 
outer wall, the boundary layer thickens to accumulate 
more low energy fluid. So the angle differs more than 
that near hub. While in other main stream regions, pas- 
sage vortices are influenced by upstream wake evidently, 
especially upper passage vortex, as the black dashed line 
region shows. 

How the flow field in nozzle blade passage is affected 
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is more interesting for designers. The detail analysis on 
shearing stress distribution on blade surfaces will be pre- 
sented. 

In Figure 7, the static pressure (SP) contour on LPT 
nozzle blade suction surface is presented with the limit- 
ing streamlines exhibited as well. Because the blade pro- 
file of all the ducts remains unmodified, the inlet flow 
angle in main stream is affected little relatively. With 
decreasing of ITD length, the boundary layer on blade 
surface shows similar information according to limiting 
streamlines distribution: radial migration towards mid- 
height passage near trailing edge remains existed. And SP 
contours almost keep unchanged: the low SP (blue) area 
locates near hub in the vicinity of throat and high SP (red) 
area sticks to leading edge. 
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Fig. 7 SP on blade suction surface 


Scrutinizing SP distribution, individual feature could 
be observed in different model to imply that the flow 
within LPT passage experiences individual accelerat- 
ing-decelerating process. In shorter ITD models, the high 
pressure area near hub decreases while it is larger near 
casing. In addition to that, low pressure area near hub 
after throat enlarges. So the radial pressure gradient is 
larger from outer wall to inner wall. On the blade surfac- 
es of ducts except L-2.2, two separation lines formed 
under the interaction between radial pressure gradient 
and passage pressure gradient could be observed. Such 
phenomenon means that the low speed fluid in boundary 
layer near walls will migrate towards mid-span height. 
However, in the passage of L-2.2, given the extremely 
large yaw angle at LPT inlet near outer casing, the sepa- 
ration line on the upper blade surface disappears or be- 
comes very weak. And the limiting streamlines are 
twisted at about 80% height, which means the upper 
passage vortex is disturbed severely by radial pressure 
gradient. 

Comparing with the flow field near suction surface, SP 


distribution on blade pressure surface is influenced more 
severely as shown in Figure 8, especially near LPT inlet. 
In shorter ducts, on one hand, the peak SP is obviously 
lower, and the continuous local high SP area in duct 
L-3.4 is interrupted to form two isolated high SP regions 
in duct L-2.2. On the other hand, as what has been men- 
tioned the inlet yaw angle below 80% passage height is 
more axial to lead larger negative attack angle. Conse- 
quently in Figure 8, a little area of separation (red dashed 
line) can be seen near leading edge in duct L-2.2 accord- 
ing to the limi- ting streamlines. What’s more, above 80% 
height, the trend of angle is different with more circum- 
ferential inlet flow in shorter ducts. As the result, the 
flow is more uneven radially to change local SP distribu- 
tion on blade surface with the decrease of duct length, 
and then the boundary layer migrates more to form a 
clearer separation line at about 85% blade height except 
L 34. 

“Clocking Effect” is an important factor for turbine 
performance, which indicates the difference of relative 
circumferential position between neighboring nozzle 
stages that could influence the flow field as well as per- 
formance of downstream turbines. Much research °'*7! 
has been done to investigate this effect. Especially 
Schennach et al.!°*! discussed the effect of upstream HPT 
nozzle wake on downstream LPT nozzle after an ITD. 
Ten models were compared to get the best relative cir- 
cumferential position between HPT nozzle blade and 
LPT nozzle blade. In this paper, similar effect is found, 
which could influence the performance of turbine as well. 
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Fig. 8 SP on blade pressure surface 


In Figure 9, Ma contours at mid-span are presented. 
The relative axial position between swirl blade and noz- 
zle blade is different as the length of duct is shortened. 
Although the circumferential positions of blades, as well 
as inlet swirl angle remain unchanged, the virtual path for 
flow is different, including upstream wake. Under vari- 
ous local pressure gradients, dissipation rate differs to 
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form individual strength and range wake (low energy) 
region in the diffuser duct. In addition to that, the relative 
circumferential position of wake in different models 
would not be definitely the same. In Figure 9, the varia- 
tions of wake circumferential position at nozzle leading 
edge caused by shortening of ITD follow evident period- 
ic feature, as the blue dashed line shows: the wake in 
L_3.4 locates just in front of nozzle blade leading edge, 
and closer to suction side; in L_3 the wake appears in the 
blade passage and closer to blade pressure side; however, 
it remains in the blade passage while locates closer to 
neighboring blade’s suction surface; finally, the wake 
migrates back near nozzle blade’s leading edge again. 
Although this periodic change of wake position is de- 
scribed according to its distribution at 50% blade height, 
it actually applies to the whole passage all the same. 
Such phenomenon is similar to modifying the relative 
azimuthal positions between neighboring stator blades, 
which would be called “Clocking-like Effect” in this pa- 
per. Such effect could change nozzle blade performance 
surely. In aero-engine, HPT rotor wake sweeps LPT noz- 
zle periodically. However, the “Clocking effect” of HPT 
blade wake should not be ignored according to literature 
B4 even through the high speed rotor blades. So the con- 
clusions reported in other papers about “the best relative 
position of nearby stator blades” should be in the same 
ITD length and inlet swirl angle, while that advice would 
be modified when ITD changes. 

Under the influence of “Clocking-like Effect”, stagna- 
tion and low velocity areas near nozzle blade leading 
edges are different. From Figure 9, low Ma area of L-3 is 
the severest, with the largest area and lowest speed. 

In addition to that, streamlines in the passage are pre- 
sented as well. Although the “Clocking-like effect” is 
obvious within the models investigated in this study, flow 
field in nozzle blade passage is well, without separation 
or back flow existing. According to Figure 8, the separa- 
tion area locates near hub in L_2.2, while flow is smooth 
in other regions. However, Ma isolines in Figure 9 exhi- 
bit individual features, under the influence of “Clock- 
ing-like Effect”. 

In order to estimate the influence of duct length on the 
loss, total pressure loss coefficient Cpt is defined: 

P-P. 
Cot = P, ti 

where Pt is the mass averaged total pressure at differ- 
ent cut planes, Pti is the mass averaged total pressure at 
duct inlet of individual duct. The larger Cpt represents 
severer pressure loss. 

Figure 10 compares Cpt at LPT nozzle outlet of dif- 
ferent models. Apparently, pressure loss caused by upper 
passage vortex (region A in Figure 10) is larger in shorter 
models. In addition to that, it breaks away from blade 
wake gradually under larger passage pressure gradient. 
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Obviously, as what has been discussed, the boundary 
layer on outer casing in shorter ducts is thicker as seen in 
Figure 5. So the leading edge horse shoe vortex near out- 
er wall would be larger. Finally, the up passage vortex, 
which is formed by pressure side branch and adjacent 
blade’s suction side branch of horse shoe vortex, would 
be larger to cause more pressure loss definitely. What 
should be analyzed again is the flow yaw angle distribu- 
tion at nozzle outlet in Figure 6. Under the effect of indi- 


Fig.9 Maat 50% passage height 
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Fig. 10 Cpt at LPT nozzle outlet 
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vidual up passage vortex, flow angle at 70%-80% pas- 
sage height differs, as the black dashed line region shows. 

Region C represents upstream wake, while being dif- 
ferent from each other in various models. In L_3.4, it is 
skewed to across almost the passage from nozzle blade 
pressure side corner to suction side neighboring blade at 
mid-height passage; in L_3 and L_2.6, it is closer to hub 
radially, even forms induce vortex with the boundary 
layer near inner wall to produce more pressure loss, 
which means wake moves along anti-clockwise direction 
in shorter models; so in L 2.2, it comes again to 
mid-height passage. This trend verifies above “Clock- 
ing-like effect” again. 

An interesting phenomenon appears for total pressure 
core region B, which is caused by lower passage vortex 
formed by leading edge horse shoe vortices near hub. 
Evidently this region presents individual feature as well. 
However, the transformation is different from that of 
region A, which is larger in shorter models. Actually it is 
the smallest in L_3, while largest in L_2.2 with separa- 
tion. The boundary layer on ITD hub is thinner according 
to Figure 5. Correspondingly the horse shoe vortices 
would be weaker in more aggressive ITD models. Such 


result is led under the influence of “Clocking-like Effect”. 


Given upstream swirl wake at nozzle outlet, it is similar 
to each other between L_3.4 and L_2.2, while similar to 
each other between L_3 and L_2.6. Especially in L_2.2, 
even a part of wake is entrained into the lower passage 
vortex to increase pressure loss evidently. 

Based on Cpt, the pressure loss coefficient of LPT 
nozzle blade is defined: 

AC = Cpto ~ Corti 

where Cpto is the C,, at LPT nozzle outlet, and Cpr is 
the Cy at inlet. 

In Figure 11, the total pressure loss coefficients of ITD 
passage (Blue-line-ITD), nozzle blade (Green line-Blade), 
as well as the total loss of ITD and nozzle blade (Red 
line-Total), are compared. 

The trend of total loss of ITD passage is almost the 
same with nozzle blade’s. And pressure loss within ITD 
duct is small. This is an interesting result that more ag- 
gressive duct would not produce more pressure loss defi- 
nitely when ITD shortens. 

In addition to that, blade passage and total pressure 
loss in shorter (more aggressive ITD) model is not larger, 
for example L-3 VS L-2.8 or L-2.6 VS L-2.4, which is 
believed affected by Clocking-like Effect. So it is rea- 
sonable to design a more aggressive ITD with larger 
radial offset and shorter length without too much increase 
of pressure loss in the LPT nozzle. The relative azimuthal 
position of LPT nozzle blade should be modified corres- 
pondingly to improve turbine’s performance further. 
However, when separation occurs (L-2.2, L-2 and L-1.8 
in this paper) the pressure loss is larger definitely. 
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Fig. 11 Total Pressure Loss Coefficient 
Conclusion 


Numerical method is applied to investigate the influ- 
ence of ITD length on LPT nozzle flow field and perfor- 
mance. Nine models are simulated. However, separation 
occurs in super-aggressive ITD models shorter than 
L_2.2, which warns the limitation of length for designers. 
Four typical models are chosen to be analyzed in detail. 
According to the discussion, several conclusions and 
advices are contributed in this paper. 

The radial angle distributions at inlet and outlet of 
nozzle of different models present the same trend to in- 
dicate good news for designers. However, at inlet the 
flow field near walls is affected more severely. The 
boundary layer at hub is thinner, while it is thicker at 
casing in shorter models. While at outlet the angle of 
main flow changes more, which implies that passage 
vortices are influenced more evidently. The shearing 
stress development on nozzle suction surface verifies 
such conclusion. Finally, the pressure loss distribution at 
nozzle outlet is affected correspondingly. 

An interesting phenomenon is “Clocking-like Effect” 
caused by the change of ITD length, which will influence 
the attack angle at leading edge of nozzle blade, even 
lead separation at pressure surface. As the result, pressure 
loss of blade passage fluctuates periodically according to 
such “Clocking-like Effect”. So it is essential to eliminate 
this effect as improving ITD through shortening length. 

In this paper, downstream nozzle flow field is found to 
be sensitive to ITD length. And separation is unavoidable 
as the ITD is too aggressive. More work about flow con- 
trol to remove separation in such super-aggressive mod- 
els is being done. 
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